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We revisit the problem of full counting statistics of particles on a segment of a one-dimensional 
gas of free fermions. Using a combination of analytical and numerical methods, we study the 
crossover between the counting of discrete particles and of the continuous particle density as a 
function of smoothing in the counting procedure. In the discrete-particle limit, the result is given 
by the Fisher-Hartwig expansion for Toeplitz determinants, while in the continuous limit we recover 
the bosonization results. This example of full counting statistics with smoothing is also related to 
orthogonality catastrophe, Fermi-edge singularity and non-equilibrium bosonization. 


Introduction .— Full counting statistics (FCS), in 
theoretical-condensed-matter context, refers to a class of 
problems involving the probability distribution of a quan¬ 
tum observable (usually the number of electrons found 
in a certain region of space or transported through the 
system over a certain time) with a particular focus on 
quantum behavior. Examples of FCS problems are the 
anti-bunching of electrons in a one-dimensional conduc¬ 
tor due to their fermionic statistics mm and a single¬ 
electron emitter, first proposed theoretically (3J and re¬ 
cently realized experimentally [4]. 

The simplest FCS problems assume non-interacting 
fermions (electrons), so that the resulting counting statis¬ 
tics can be expressed in terms of a determinant taking 
into account antisymmetrization of the relevant multi¬ 
particle processes |T|. An alternative approach is based 
on the bosonization technique (using the equivalence 
between bosons and fermions in one dimension) [6|. 
Bosonization methods for FCS can be extended to in¬ 
clude interactions, but they usually do not fully take into 
account the discreteness of particles HEHg. 

It is therefore important to understand limitations of 
the bosonization approach to FCS and its connection to 
the exact calculation with discrete fermions. In this Let¬ 
ter, we address this problem by studying FCS in a one¬ 
dimensional free-fermion model, where an uncertainty 
(smoothing) in the counting procedure is introduced, so 
that the particle number is no longer quantized. In this 
model, we can study in full detail a crossover between 
the discrete-particle and bosonization results as the un¬ 
certainty increases. While the details of this crossover 
depend on the profile of the introduced uncertainty, the 
qualitative description of the crossover is found to be uni¬ 
versal: FCS respects the discreteness of particles if the 
uncertainty region is much narrower than the (average) 
inter-particle distance, but crosses over to the bosoniza¬ 
tion result when the uncertainty region is much wider 
than the inter-particle distance. 

Mathematically, this problem amounts to studying 



0 10 x 


FIG. 1: The weight functions wl,o.{x) for L = 10 and a = 1 
for the three models of smoothing. 


an evolution of the asymptotic behavior of a Toeplitz 
(or, more precisely, Wiener-Hopf) determinant with a 
Fisher- Hartwig singularity mm as the singularity is 
smoothed in a certain way. In the case of a sharp singu¬ 
larity, the corresponding determinant is given by a double 
asymptotic series of a Fisher-Hartwig type fully respect¬ 
ing the particle discreteness HBH2Q1. As the singularity is 
smoothed, the secondary branches of this expansion get 
suppressed, and the remaining leading branch reproduces 
the bosonization result. 

Model .— We consider the one-dimensional model of 
spinless free fermions on a line (both coordinate and mo¬ 
mentum are continuous) at zero temperature. The sys¬ 
tem is in its ground state, which is the Slater determinant 
of plane waves characterized by the Fermi wave vector 
kF' the states with the wave vectors smaller than fep 
(in absolute value) are occupied, and the states with the 
wave vectors larger than kF are empty. We introduce the 
FCS generating function as mm 

XL,a{n) = (exp(27T inQ L , a )) , 


(1) 
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where the average is taken over the ground state. The 
particle-counting operator Ql,o. is defined as 

/ OO 

'i> i {x)'^(x)w Lj a{x)dx, (2) 

-OO 

where ^(x) and d'(x) are the fermionic creation and 
annihilation operators and w F , a (x) is the weight function 
for the particle counting. The weight function depends 
on the two parameters: the length of the interval L and 
the smoothing length scale a, Fig. [l] We always assume 
L > a. 

In the original formulation of the FCS problem, the 
smoothing is absent (a = 0), and the weight function is 
the characteristic function of a line segment, 


w L , o(x) = 6(x)0(L - x) 


0 x < 0 or x > L, 

(3) 

1 0 < x < L. w 


In this case, Ql.o. takes integer values, and the gener¬ 
ating function (jT]) is periodic in the counting variable, 
Xl, o(k + 1) = Xl, o(«). It can be expressed as a Toeplitz 
determinant (see Appendix), which has been subject to 
extensive studies HM7]. The periodicity of xl, o(k) is re¬ 
flected in a Fisher-Hartwig-type asymptotic series (con¬ 
jectured and numerically verified to a very high order) 

EEH2U: 

OO 

XL,o(k)= Xl,o(k+j) , (4a) 

j=- OO 


XL, o(«) = exp 


2'kikNl 


2k 2 ln(2nN L ) 


+ C(n) + J2fn(K) ( iN L )~ n 

71=1 


(4b) 


The coefficients in this double expansion can be explic¬ 
itly calculated: C(k) is expressed in terms of Barnes G 
functions as C(k) = 21n |G(1 + re)G(l — k)|, and /„(k) 
are polynomials in n, which can be computed iteratively, 
order by order [2D] (the same coefficients are denoted 
P nn { k) in [21]). xl, o(k) depends only on one parameter: 
the average particle number on the segment, 


N l = k F L/n . (5) 

On the other hand, a bosonization approach proposed 
for the generating function xl,o,(k) [2j leads to 

XL°a( K ) = ex P ^ttzkNl-2k 2 [\n(N L /N a ) + const]) , (6) 

where we have defined the average number of particles in 
the smoothing region 



FIG. 2: (a)—(c): Coefficients Co(k, a) and Ci(k, a) as a func¬ 
tion of k for several values of N a , as computed from numeri¬ 
cal fits. Co(k, a) are shown by solid lines, Ci(k,ci) by dashed 
lines. Different symbols correspond to different values of N a 
as shown in the legends. In panel (c), the dash-dotted line is 
the bosonization result (12 I at N a = 0.125. (d): Coefficient 
C i(/t=l, a) as a function of N a for the three models. The plot 
for Lo model is the exact analytical result (201. The dotted 
line is the leading asymptotics (191 for RC model. 


and the constant at the logarithm depends on the details 
of the smoothing. The bosonization expression ([b]) corre¬ 
sponds to the leading terms in the leading branch (j = 0) 
of the double expansion (§• We expect (and confirm it 
the present work) that this result applies provided the 
smoothing a is sufficiently large. It is not periodic in k, 
which is consistent with a continuous spectrum of the 
operator Ql,o. at a finite a. 

The goal of the present study is to examine in detail 
the crossover from the discrete particle counting given 
by the expansion Q to the bosonization result © as 
a function of the smoothing parameter a. Some of the 
results are derived analytically, while a more detailed pic¬ 
ture of the crossover is obtained numerically. Smoothing 
is introduced as a convolution 
r°° dx' 

w L ,a{x) = — 6(x)8(L - x) g(x'/a ), (8) 

J — OO ® 

where g(x) is a smoothing function for the scaled coor¬ 
dinate x = x/a. We always assume the normalization 
f dx g{x) = 1. To be specific, we consider three models 
of smoothing corresponding to the three choices of the 
smoothing function: RC (this type of smoothing may be 
relevant for particle-counting in time, where the counter 
has a RC-type response function), Lorentzian (Lo) and 
Gaussian (Ga): 

(0(x) exp(— x) RC, 
g(x) = < [tt(x 2 +1)]" 1 Lo, (9) 

[exp(— x 2 / 2 )/v / 27 t Ga 


N a = k F o J /'K 


( 7 ) 
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(see Fig. [I]). In any of these models, the FCS generating 
function x_L j0 (re) depends on L and a via the two dimen¬ 
sionless parameters Nl and N a . 

Results .— We conjecture (and support this conjec¬ 
ture by analytical and numerical results) that, at a finite 
a, the FCS generating function admits the same struc¬ 
ture of the expansion Q as at a = 0, except that the 
expansion coefficients are not periodic in the branch in¬ 
dex j and depend on the smoothing parameter a. If we 
limit the range of re to the interval [0,1], we can write the 
leading terms of the expansion as 


XlA k ) ~ ex P 


2n ikNl 


2k 2 In Nl + Co (re, a) 


where 


T = 7 + lim 
£—>-0 



dk 

~r9k g-k + lne 


(13) 


in the case of a general smoothing model ([8|. Here 7 = 
0.5772... is the Euler-Mascheroni constant and g k are 
the Fourier components of the smoothing function, 


/ OO 

dxe~ lkx g{x)- (14) 

-OO 

For the three models considered in our paper, a calcula¬ 
tion gives 


+ exp 


27 n(re 


1)N l - 2(re - l) 2 lnN L + Ci(k, a) . 

( 10 ) 


At a = 0, in agreement with the full expansion the 
coefficients Co(re,a) and Ci(re,a) are given by 


C 0 (re,0) =<7 i(1-k,0) 

= 2 In |G(1 + re)G( 1 - re)| - 2 k 2 ln(27r). (11) 

We expect that, as a increases, Co(re, 0) also increases 
to reproduce the bosonization result ([6]), while Gi(re, 0) 
decreases to suppress the second term in the expansion 

This behavior is indeed confirmed by numerical stud¬ 
ies. We report the details of our numerics in the Ap¬ 
pendix, an here we only plot in Fig.[2£i-c the results of the 
numerical fits for the coefficients Go (re, a) and Gi(re,a) 
for the three models of the smoothing. Figure [2]l shows 
more detailed results for the coefficient Ci(l,a) repre¬ 
senting the suppression of the secondary Fisher-Hartwig 
branch at re = 1. This coefficient was, in fact, com¬ 
puted as a doubled contribution from a single “phase 
slip” (at the beginning or at the end of the interval) in 
the weight function wlA x ) ( n °t e that at re — 1 the val¬ 
ues WL, a {x) = 0 and Wl,u{ x ) = 1 are equivalent, see 
Appendix). An agreement between the results computed 
with this method (at re = 1) and the fitting procedure for 
the full function \L,a { K ) serves as an independent check 
of our numerical scheme. 

As expected, while the actual values of the coefficients 
Go(re,a) and C i(re, a) depend on the chosen smoothing 
model, qualitatively the crossover between the discrete 
particle counting and a continuous bosonization descrip¬ 
tion occurs at N a ~ 1 in all the three models. 

Furthermore, some of the numerical results presented 
above can be verified by analytical means. First, the con¬ 
stant in the bosonization formula © can be computed 
from the corresponding Toeplitz determinant using the 
strong Szego theorem [22; [ 23 ] . This results in the asymp¬ 
totic behavior for Go(re,a): 

Co (re, a —> 00 ) = 2re 2 (ln N a — T), (12) 


Trc = 7 1 


T Lo = — In 2 , 


T Ga = 7 /2. (15) 


These analytical results perfectly agree with our numer¬ 
ical fits: for the RC and Lorentzian smoothing, the com¬ 
puted values of Co(re, a) for N a > 0.25 and N a > 0.125, 
respectively, are indistinguishable in the plots of Fig. §'-b 
from the bosonization asymptotics (12), (15). For the 


Gaussian smoothing, the difference between Co (re, a) and 
the bosonization asymptotics is visible for N a = 0.125, 
but not for N a > 0.25 (Fig. [2J3). 

Second, the asymptotic behavior of Ci(l, a) as a —> 00 
can also be calculated analytically using an asymptotic 
formula for the Toeplitz determinant with a non-zero 
winding number [24]. A calculation results in 


-Ci(l,a 


00) 

dz 


+ Re In / — exp 

27T 


- T 

2iriN a z — 


dk g k 


1 — e 


ikz " 


|fc| 


where T is given by (131 and S is defined as 


= 2y + lim 
£—>-0 


/ °° rib 

T (g k + g _ k ) + 2 lne 


(16) 


(17) 


An explicit calculation gives 
Src = 27 , S Lo = 0 , 


lGa = 7 + ln2 - ( 18 ) 


Remarkably, the asymptotics (16), while not exactly co¬ 
inciding with Ci(l,a), gives a very good approximation 
in the whole range of the values of a (in the plot in Fig. ©1 
it would be indistinguishable from the exact values, see 
details in Appendix). 

The analytic approximation (16) also allows us to 
extract the leading asymptotic behavior. For the RC 
model, we find 


Ci(l,a —»• oo)r C ~ 27 — 41n(27rlV a ), (19) 

i.e., the secondary branch of y(l, a) decays as N~ 4 . For 
the Lorenzian smoothing, the Toeplitz determinant for 
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a single “phase slip” w a (x) = ^ hi(jr^) may be com¬ 
puted exactly by using the analyticity of the function 
exp[27rw Q (a;)] in one of the half-planes. An argument in 
the spirit of EU25] then leads to the exact result: 


C , i(l,o) Lo = - 47 rJV„. 


( 20 ) 


Finally, for the Gaussian smoothing, the asymptotic be¬ 
havior is 

Ci(l, a -» oo) Ga ~ -47riV a v / ln(27rAT2) _ ( 21 ) 

However, this expression by itself [unlike the integral for¬ 
mula (16)] does not provide a good approximation for 


Ci(l,a), since the next-order corrections are not 0 ( 1 ). 

Remarkably, even though the coefficient at the sec¬ 
ondary branch is reduced with increasing N ai the depen¬ 
dence on N l does not change. In particular, for k > 1/2, 
the bosonization (first) term in ( 10 ) always decreases 
with Np faster than the non-bosonization (second) one. 
This implies that the non-bosonization corrections aris¬ 
ing from the discreteness of particles should be visible, 
in this range of k. at sufficiently large L (see Appendix). 
The corresponding values of L depend on the particular 
choice of smoothing and can be deduced from comparing 


the two terms in (10) using the asymptotics at large N a 


derived above. A practical observation of this effect may 


be limited at large N a if both terms in (10) become too 
small. 

Discussion In this Letter, we have discussed the 
crossover from discrete to continuous FCS in the model 
of free one-dimensional fermions, as a function of smooth¬ 
ing. This model may serve as an illustration of a relation 
between discrete and continuous descriptions in a wide 
spectrum of similar problems, including FCS with tempo¬ 
ral measurements, Fermi-edge singularity, orthogonality 
catastrophe, and non-equilibrium bosonization. 

By FCS with temporal measurements we understand a 
FCS setup where the counting is performed over a certain 
time interval, e.g., by opening and closing an electric con¬ 
tact or by applying a time-dependent voltage pulse mm- 
This formulation of FCS used in the original FCS papers 
differs from our approach where the measurement is ex¬ 
tended in space instead of time. This difference implies 
the necessity of regularizing the contribution of the Fermi 
sea in temporal FCS: a large current of left-movers is 
compensated by a similarly large current of right-movers, 
so the average charge transfer is determined by the vicin¬ 
ity of the Fermi level, while the quantum noise depends 
on the ultraviolet cutoff defined by the Fermi energy Ep 
(arising from the bottom of the Fermi sea). In our spatial 
formulation, there is no cutoff introduced by the bottom 
of the Fermi sea, and both the total number of parti¬ 
cles and its fluctuations depend on the same scale kp. 
As a consequence of this difference in ultraviolet cutoffs, 
our results cannot be literally translated to the case of 
temporal FCS. However, we expect that the main quali¬ 
tative conclusion will hold also in the temporal case: the 


generating function x( K ) loses periodicity in k when the 
smoothing time scale exceeds the typical time between 
particles (i.e., s~j ? 1 ). 

Determinants similar to the generating function x( K ) 
also appear in problems related to the orthogonality 
catastrophe and to the Fermi-edge singularity (FES). In 
those cases, the counterparts of the secondary branches 
of x( K ) are secondary singularities (cusps or peaks) in 
the frequency-dependent response function (the closed- 
loop contribution in the FES context). Such secondary 
cusps and peaks were studied, e.g., in the recent work [26] 
(and, in the case of a bound-state contribution, earlier in 
[27]). According to our predictions, such cusps and peaks 
should be most visible in case of instant switching of the 
scattering potential, but get suppressed if the switching 
time of the scattering potential exceeds ef . 

We should also remark that most studies of FES in¬ 
volve an artificial regularization of the Fermi-sea con¬ 
tribution at energy scales smaller than Ep, so that sec¬ 
ondary singularities in the response function are ne¬ 
glected [25H50] . Such secondary singularities (separated 
from the main peak by ep) are probably not experimen¬ 
tally relevant in physical metals, where Ep is a large en¬ 
ergy scale, but may be of interest in other models with 
FES physics where Ep is, for some reason, small (see, e.g., 
Ref. [H] for an example from spin-liquid theory). 

Finally, we also mention recent works on non¬ 
equilibrium bosonization where determinants similar to 
ours appear (again, in the temporal form) [3]. Similarly 
to the FES problem discussed above, these works assume 
a regularization of the Fermi-sea contribution, which is 
equivalent to neglecting the secondary branches of the 
Fisher- Hartwig expansion related to the bottom of the 
Fermi sea. 
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APPENDIX 

A. Toeplitz determinant for x{k )■— The FCS gener¬ 
ating function (JT|) can be written as a trace in the multi- 
particle space and then re-expressed as a determinant in 
the single-particle space mmm- 

XL,ai K ) = det [(! - n F) + n F e 2mKW ] , (22) 

where np is the Fermi occupation number and w is the 
operator of multiplication by WL, a {x). The operator np 
is diagonal in the momentum space, while w is diagonal in 
the coordinate space. In this paper, we consider the zero- 
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FIG. 3: (a) the values of Ci(re=l,a) for RC and Ga mod¬ 
els (the same as in Fig. gl) as calculated numeri call y, to¬ 
gether with the analytic approximations given by ( |16| ). (b) 
the difference between the approximate and exact values of 
Ci(k=1,o) (upper curve: Ga model, lower curve: RC model). 


FIG. 4: (a) In |x'L,a(«)| as a function of k in the RC model 

at Nl = 10 for selected values of N a . Periodicity in k is lost 
with increasing N a . (b) ln|xz,,a(re)| as a function of Nl in 
the RC model at k = 1 and N a = 1.5. For large Nl, the 
second term in (101 dominates. 


temperature case, so hf is a projector onto the occupied 
states (below fc_p). 

In the limit of no smoothing (a = 0), w is also a pro¬ 
jector (in the coordinate space), and the determinant 
is symmetric with respect to exchanging coordinate and 
momentum: 


Xl, o(«) = det [l - (1 - e 2mK )n F w\ 


(23) 


This operator is of Toeplitz (or Wiener-Hopf) form in 
both coordinate and momentum representations. 


After introducing smoothing, the operator in (22) is no 


longer Toeplitz in the coordinate representation, but re¬ 
mains Toeplitz in the momentum representation (at zero 
temperature), 

B. Toeplitz determinant for Ci(re=l ,a ).— At re = 1, 


the term with C\ (re, a) in the decomposition (10) does not 
oscillate with L. It therefore arises from the end points 
of the segment of measurement (x = 0 and x = L). This 
allows us to calculate Ci(re=l, a) directly from a single- 
step contribution: 


-Ci(re=l,a) = Re In det [(1 — njr) + n^e 




(24) 


where w is the operator of multiplication by the single- 
step counterpart of WL,a{x ): 


.(*) = f 

J — C 


dx' 


9{x) g{x'/a). 


(25) 


Note that while w a (x) has different limits at x —> ±oo, 
the exponent exp(27mD) tends to 1 in both limits, and 
the determinant ( f2lj ) is well defined. 

C. Details of the numerical calculation .— For nu¬ 
merical calculations of Xl,u{k ) and Ci(re=l,a), we ap¬ 
proximate the integral operators (22) and (24) by finite¬ 
dimensional matrices by discretizing the momentum 
space (which is equivalent to considering the system on 
a circle). In the momentum space, these matrices are of 
Toeplitz form. The determinants of those matrices are 


calculated numerically for several matrix sizes and then 
extrapolated to the infinite matrix size in order to obtain 
XL,a(n) and Ci(re=l,a), respectively. 

For numerical calculations of the determinants, we use 
the LAPACK library [ 33 ], and Fourier transformations 
necessary for computing the matrix elements were done 
with the help of the GNU scientific library [ 31 ]. In our 
calculations, determinants of matrices of the linear size 
up to 600 were calculated, which allowed us to achieve 
good precision for the considered ranges of parameters. 

The values of XL,a{x ) were then fitted to the expan¬ 
sion ([To]) , and the coefficients Co (re, a) and Ci(re, a) were 


extracted (the fit also included terms up to L~ 2 in the 
exponents). As an additional check of this fitting pro¬ 
cedure, we have verified that the values of Ci(re=l,a) 
obtained from these fits agree with those calculated in¬ 
dependently using the one-step method described above. 

D. Analytic approximation for Ci(re=l, a—>oo ).— 


The approximation (161, while derived under the assump¬ 


tion N a 1, is remarkably accurate even at small N a . In 


the Lo model, where Ci(re=l,a) is given exactly by (20), 


the approximation (16) reproduces the same exact result 


for all N a . In the RC and Ga models, the difference be¬ 
tween (16) and Ci(re=l,a) tends to zero as a —» oo, but 
remains finite at finite a. However, this difference is nu¬ 
merically small for these two models, even at a = 0. We 
plot this difference as a function of N a in Fig. [ 3 ] 

E. Numerical illustration of the crossover .— We il¬ 
lustrate the crossover described in the main body of the 
paper by two plots. In Fig. |4^, we show how the peri¬ 
odicity of Xi,a( K ) i n K is gradually lost as N a increases. 
Figure |4 ]d illustrates how the second term in (10) domi¬ 
nates at large L (in RC model, for re = 1, this happens 
at kpa 3 ). 
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